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Many  studies  have  been  made  concerning  the  use  of 
i,  iriant— imbedding  principles  to  obtain  integrcr-differential 
uations  describing  transport  processes.  In  this  Memorandum 
M  author  formulates  particle-transport  processes  in  terms 
m  abstract  mathematical  system.  Thia  study  evolved  from 
research  in  neutron-transport  theory  sponsored  by  the 
Advanced  Research  Projects  Agency. 

This  version  was  rewritten  in  order  to  achieve  greater 
clarity,  and  it  diff  rs  from  the  original  version  by  a 
change  in  notation  and  emphasis. 


A  complete  mathematical  definition  of  an  abstract 
linear  transport  process  is  given  in  terms  of  a  new 
axiomatic  system.  After  several  preliminary  theorems  are 
proven,  the  basic  algebraic  equations  relacing  the 
transport  operators"  are  derived,  and  the  semigroup 
properties  of  these  equations  are  indicated.  These 
algebraic  equations  are  then  used  to  derive  the  standard 
differential  equations  describing  a  transport  process. 

The  generalized  irvariant-imbedding  equations  are 
used  to  describe  an  energy-dependent,  neutron-transport 
process. 
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A  NEW  APPROACH  TO  TRANSPORT  THEORY 


1.  INTRODUCTION 

The  methods  of  invariant  imbedding  [1]  have  been 
utilized  by  Bellman,  Kalaba,  and  Wing  [1,2,3]  to  analyze 
one-dimensional  neutron— transport  processes  arising  from 
the  chain  reactions  of  nuclear— fission  processes.  The 
invariant  imbedding  technique  allows  one  to  replace  the 
classical  equations  of  transfer,  which  are  constrained  by 
two— point  boundary  conditions,  with  a  system  of  differential 
equation;?  with  initial-value  constraints.  The  relation 
between  the  invariant-imbedding  approach  to  neutron  transport 
theory  and  the  circuit— theory  formalism  of  transmission 
lines  [4]  was  discovered  by  Redheffer  [4,5,6],  who  investi¬ 
gated  the  algebraic  structure  of  transfer  processes. 
(Extensive  references  on  neutron  transport  and  transmission 
lines  can  be  found  in  [3]  and  [5].)  In  this  Memorandum,  we 
propose  a  mathematical  formulation  of  one— dimens ional 
transport  processes  which  yields  the  algebraic  description 
of  the  process  in  terms  of  Redheffer1 8  "scattering  matrix" 

[5]  as  well  as  the  differential  equations  arising  from 
Bellman's  invariant— imbedding  approach  to  neutron  transport. 

A  result  of  our  mathematical  description  is  a  demonstration 
of  the  nonsingularity  of  the  "transmission  operators" 

(See  c.  2)  which  is  independent  of  any  assumption  of 
differentialability  of  the  process.  We  also  apply  our 
results  to  obtain  integro-differential  equations  which 


completely  describe  one-dimensional  energy— dependent  neutron 
transport  processes. 

We  shall  let  the  one-dimensional  neutron  transport 
process  serve  as  the  motivating  example  for  our  mathematical 
system.  The  one-dimensional  transport  process  consists  of 
a  finite  real  line  segment,  referred  to  as  the  rod,  which 
represents  the  interacting  medium  through  which  the  neutrons 
tnvel  (see  Fig.  1).  The  neutrons,  considered  as  abstract 
point-particles,  travel  through  the  rod,  moving  either  to 
the  right  or  to  the  left,  and  interact  with  the  medium, 
occasionally  producing  additional  neutrons  by  nuclear 
fissions.  For  each  point  y  of  the  rod,  we  let  u(y)  and 
v(y)  represent  the  expected  flux  of  neutrons  passing  through 
y  to  the  right  and  to  the  left,  respectively  (see  Fig.  1  ) . 

— -  v<» 

- -  u(yj 

I - 1 - 1 

0  y  x 

Fig.  1-Rod  Extending  from  0  to  x 

To  keep  our  discussion  general,  we  shall  assume  that  the 
fluxes  u(y)  and  v(y)  are  elements  of  a  Banach  space.  The 
transport  process  can  then  be  described  in  terms  of  linear 
differential  equations  of  the  form 

(i.i)  u<y>  •  bi<y)u(y)  +  *(y)v(y) , 

fly  v(y)  »  -c(y)u(y)  -  b2(y)v(y)  , 
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where  u(y)  and  v(y)  are  the  fluxes  to  the  right  and  left 
at  the  point  y,  and  a(y)  ,b^(y)  ^(y) ,  and  c(y)  are  linear 
operators,  defined  for  each  point  y  in  the  rod  [0,x]. 

The  formulation  of  transport  processes  in  terms  of 
equations  (1.1)  has  a  major  drawback.  To  study  practical 
problems  in  which  incident  neutron  beams  initiate  a  trans¬ 
port  process  in  the  rod  [0,x],  we  must  solve  equation  (1.1) 
subject  to  the  boundary  conditions 

u(°)  -  uin,  v(x)  -  vin  , 

where  u^n  and  v^n  are  the  incident  fluxes  to  the  right  at 
0  and  to  the  left  at  x,  respectively.  (In  many  applications 
we  pick  uin  or  v^n  to  be  0.)  Thus  we  must  solve  equations 
of  the  form  (1.1) — which  may  be  2n  simultaneous  ordinary 
differential  equations  [3]  or  integro-differential  equations 
[2] — with  two-point  boundary— value  conditions.  This 
formulation  is  often  unsatisfactory  when  numerical  answers 
are  desired,  because  of  the  difficulty  of  numerically 
solving  two-point  boundary— value  problems  of  high  dimen¬ 
sionality. 

The  technique  of  invariant  imbedding  has  been  success¬ 
fully  used  by  Bellman,  Kalaba,  and  Wing  [1,2]  to  reduce 
the  two-point  boundary— value  equations  of  neutron-transport 
processes  to  initial— value  problems.  We  now  briefly  describe 
the  invariant— imbedding  approach.  Consider  a  one-dimensional 
transport  process  on  the  rod  [0,  L]  which  is  described  by 


the  linear  operators  a(y),  b^(y),  b2<y)i  c(y)  (for  y  e  [0,L]) 
of  equation  (1.1).  For  each  x  in  the  interval  [0,L],  we 
consider  the  transport  process  taking  place  on  the  shorter 
rod  (0,  x]  which  is  described  by  the  same  linear  operators. 
We  define  the  family  of  operators  P(x),  r(x),  t(x),  and  t(x) 
as  follows: 

p(x)  f  *»  the  "reflected"  flux  leaving  the 
rod  [0,  x]  to  the  right  due  to  a 
flux  $  incident  upon  the  rod  to  the 
left, 

t(x)  ♦  ■  the  "trar emitted"  flux  leaving  the 
rod  [0,xj  to  the  left  due  to  a  flux 
ijr  incident  upon  the  rod  to  the  left, 
r(x)  cp  *  the  "reflected"  flux  leaving  the 

rod  [0,x]  to  the  left  due  to  a  flux 
cp  incident  upon  the  rod  to  the  right, 
t(x)  cp  -  the  "transmitted"  flux  leaving  the  rod 
[G,xJ  to  the  right  due  to  a  flux  cp 
incident  upon  the  rod  to  the  right. 

Thus  in  the  particle-transport  process  on  [0,x],  If  the 
incident  fluxes  are  given  by 

u(0)  -  cp,  v(x)  -  ♦, 

then  the  "output"  fluxes  (see  Fig.  2)  are  given  by 
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u(x)  ■  t(x}  Cp  +  p(x)  ♦, 
v(0)  ■  r(x)  cp  +  t(x)  i|r. 


r(x)cp  +  t(x)i|i  - - 

cp  - 

h- 

0 


- —  ♦ 

- t(x)cp  +  P(x)i|i 

H 

x 


Fig.  2 

By  considering  the  rod  [0,  x  +  A],  and  using  the  "particle 
counting"  technique  of  Bellman  et  al.  (see  [1,  2,  3]),  we 
arrive  at  the  operator  differential  equations 

(1.2)  gj  P(x)  -  a(x)  +  b1(x)p(x)  +  p(x)b2(x)  +  p(x)c(x)p(x), 

—  t(x)  "  t(x)  [b«(x)  +  c(x)P(x)], 

<L' 

—  r(x)  -  t  (x)  c  (x)  t  (x)  t 
dx 

—  t(x)  *  [b1(x)  +  P(x)c(x)]t(x), 
dx 

with  the  initial  conditions 

(1.3)  P( 0)  -  r(0)  -  0  , 
t(0)  -  t(0)  -  1  , 

where  I  is  the  identity  operator. 

The  operators  P,r,T,  and  t,  called  the  transport 
operators .  play  an  important  role  in  the  theory. 


By  using  the  transport  operators,  one  replaces  the  two- 
point  boundary-value  problems  of  classical  transport  heory 
by  the  differential  equation  \1.2)  with  the  initial  conditions 
(1.3). 

2.  THE  MATH31ATICAL  FORMULATION 

We  now  give  an  axiomatic  definition  of  transport 
processes  from  which  we  shall  mathematically  derive  our 
desired  results.  To  account  for  the  ''infinite-dimensional 
case,"  we  shall  consider  the  transport  operators  as  linear 
operators  on  a  Banach  space.  We  shall  closely  follow  the 
notation  used  by  Redheffer  [5,6]. 

Lotjfbe  a  Banach  space.,  over  a  real  or  a  complex 
field.  We  shall  refer  to  elements  oijQ  as  fluxes  in  order 
to  recall  the  neutron  transport  process  discussed  in  the 
introduction.  (When  dealing  with  transmission  lines, 
elements  of  V  would  represent  electrical  signals.)  We 
shall  also  consider  the  Banach  sp  ce  ■  X 
shall  write  elements  of  X  ^  in  the  form  (u,v) 
u,v  cX  •  The  norm  1x1  X2  is  given  by 

1 1  (u,v)  1 1  -  Null  +  I  Ml, 

where  | |uj |  is  the  given  norm  in  X  * 

Defini  ion  1.  Let  the  real  interval  [Xq,x^J  be  given. 
Define  the  space 
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2 

considered  as  a  subspace  of  the  metric  space  R  . 

Let  S  be  a  mapping  of  D  into  the  space  of  bainded  linear 
operators  onjC^’  Thus  for  every  (ypy2)  w®  associate 

with  it  a  bounded  linear  operator 


s(y!,y2):  j}  "*  X2 


We  also  write 


u 


(2.1)  S(yi,y2) 


t(ypy2)u  +  P(y1iy2)v 

r(y1,y2)u  +  T<ypy2>v 


ft(y!^y2)  p(yry2> 

r(yi>y2)  Tiy^) 


u 


for  all 


x2. 


*-  — J 

Since  S(y^,y2)  is  always  assumed  to  be  a  bounded  linear 

A 

operator  on£  ,  it  follows  that  t(y1,y2),  T(y1,y2),  r(yl,y2) , 
and  P(y^#y2)  are  bounded  linear  operators  onjC  .  We  can 


abbreviate  (2.1)  by  writing 

S  - 


rt  ? 


We  pause  to  describe  the  physical  meaning  that  we 
attach  to  the  above  definition.  Considering  the  neutron- 
transport  model  mentioned  in  the  introduction,  we  regard 
each  element  (y^y2)  of  D  as  the  rod  segment  extending 
frcm  y^  to  y2  "imbedded”  in  the  whole  rod  [x^x^].  The 
linear  operator  S(y^,y2)  represents  the  transformation 
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relating  the  "output  fluxes"  from  the  rod  [y^^J  to  the 
"incident  fluxes."  Specifically,  we  want  to  have 


(2.2) 


“(y2^  e,  S  u<y,) 

S(yl'y2>  v(y2)  ’ 


where  u(yj)  and  v(yt)  are  the  fluxes  to  the  right  and  to 
the  left,  respectively,  at  the  point  yp  i  -  1,2  (see  Fig.  3) 
The  reader  should  verify  that  t(yi,y2),  T<yi,y2),  r(ypy2), 
*nd  P(ypy2)  represent  the  transport  operators  (as  defined 
in  Sec.  1)  for  the  rod  [ypy2J.  The  operator  S(ypy2)  is 
called  the  scattering  operator  in  analogy  with  Redheffer's 
scattering  matrix  [5,6]. 


u(yx) 

v(yj) 


n(y2) 

v(y2) 


Fig.  3 


Definition  2.  For  each  (y1,y2>€D,  define  the  relation 
G*1  c  X2  x  X2  by 

y2 

yl  f 

Gy  "  |(w].*w2)  I  (u^vp  *  S(y1,y2)(u1,v2), 


where 


.'vi>  1-1.2}  . 


Each  element  of  is  thus  an  ordered  pair,  the  first 

y2 

dement  representing  the  fluxes  to  the  right  and  left* 
respectively, at  the  point  yp  and  the  second  element 
representing  the  fluxes  at  the  point  y2*  We  say  that 
the  relation  G  relates  the  fluxes  at  the  point  y0  to 

y2  * 

the  fluxes  at  the  point  yp  in  the  transport  process  under 

consideration.  Note  that  we  are  not  assuming  that  given 

the  fluxes  u(y^)  and  v(yj)  at  the  point  yp  there  would 

exist  fluxes  u(yj)  and  v(yj)  satisfying  equation  (2.2). 

In  the  course  of  our  mathematical  analysis,  we  shall  instead 

show  that  there  is  a  one— one  relationship  between  the 

fluxes  at  two  different  points,  and  equations  (2.3)  and  (2.4) 

give  an  explicit  formula  for  the  fluxes  at  one  point  in 

terms  of  the  fluxes  at  another  point.  Until  we  demonstrate 

this  one-one  correspondence,  however,  we  shall  have  to  be 

content  with  the  somewhat  clumsy  representation  in  terms  of 

yi 

the  family  of  relations  G  . 

y2 

We  now  review  briefly  pome  basic  mathematical  con¬ 
ventions  to  be  used  in  the  following  discussion.  The 
space  of  bounded  linear  operators  is  given  the  topology 
induced  by  the  usual  norm, 

llTll  -  ,8UP  I  |Tw|  I . 

|  |w|  1-1 

(This  is  called  the  uniform  operator  topology.)  A  linear 
operator  on  a  linear  space  is  said  to  be  invertible  when  it  is 
one-one  and  onto.  If  T  is  an  invertible  bounded  linear 
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operator  on  a  Banach  space,  then  the  elementary  theory  of 

—1 

Banach  spaces  [7]  tells  us  that  the  inverse  T  is  also 
bounded.  The  composition  of  two  relations  and  G2  is 
defined  by 


G2G1  “  Y  such  that  (oijYjcG^  and  (y,3>€G2} 


If  the  relations  G^  and  G2  are  functions,  then  their 
composition,  G2G^  is  the  usual  composition  of  functions. 

We  are  now  ready  to  give  a  complete  definition  of  a 
one-dimensional  transport  process. 

Definition  3.  Let  S  be  described  as  in  Definition  1. 

We  say  that  S  describes  a  transport  process  when  the 
following  three  axioms  are  satisfied: 

(i)  S(y>y)  -  I  for  all  y  e  [xq^J 
(I  is  the  identity  operator) ; 

(ii)  S  is  continuous  in  the  uniform  operator 
topology i 

(ili)  gY3  Gy\  ■  S3  £or  *0  <  yl  £  y2  *  *3  £  *1' 

where  G  is  given  by  Definition  2. 

yj 

In  the  remaining  discussion  we  shall  always  assume 
that  S  is  given  and  describes  a  transport  process. 

Axiom  (i)  is  equivalent  to  the  statement  that 
r(y,y)  -  P(y,y)  ■  0  and  t(y,y)  -  r(y,y)  -  I.  Axiom  (iii) 
is  really  a  shorthand  notation  for  the  following  two  .  ^tements: 


a. 
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(a)  If 


then 


(b)  If 


then 


fu# 

LviJ 


u. 


u. 


1 

"u2 

v2 

S(y1,y2) 

sCypya) 

s(yvy%) 

.  ,2 


u. 


2 

“i 


fui 


It. 


3' 


u. 


such  that 


u. 


s(y1yy2) 


u. 


fu. 


and 


s(y2,y3) 


and 


“31 
J  i 

LV2j 


-  S(y2,y3) 


u2 

v3 


u. 


Figure  4  serves  as  an  illustration  of  this  cardinal 
axiom.  Note  that  we  do  not  assume  that,  for  a  given 
transport  process  in  a  rod  [y^y^]*  there  exist  unique 
internal  fluxes  consistent  with  the  process  at  an  inter¬ 
mediate  point  y2>  We  shall  instead  prove  this  uniqueness 
in  the  course  of  our  mathematical  analysis. 
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Lemma_l:  3A  >  0  such  that  h  €  [0,Aj  and  (y,y  +  h)  e  D 
implies  that  t(y,y  +  h)  and  T(y,  y  +  h)  are  invertible. 
Proof.  Let  tCy^^)  -  I  -  Z(y1,y2>.  Since 


I |z(y|,  y2)l I  “  I l2(y1,y2) | I 


£  1 1  t(yx,y2) 


t(y1^y2) i I 


1  I  |S(ypy2)  -  S(y1#y2)| 

the  continuity  of  S  (in  the  uniform  operator  topology) 
implies  the  continuity  of  ||Z| |.  Since  D  is  compact,  we 
conclude  that  ||Z||  is  uniformly  continuous  on  D.  There¬ 
fore  ^  >  0  such  that  h  e  [0,6^]  and  (y,y  +  h)  e  D  implies 

that  | |Z(y,y  +  h) | |  <  1,  and  therefore  that 
t(y,  y  +  h)  -  I  -  Z(y,  y  +  h)  is  invertible  (see  [7], 
page  66).  Similarly,  3  $2  such  that  he[0,62j  and  (y,  y  +  h)  e  D 
implies  that  r(y,  y  +  h)  is  invertible.  We  then  pick 


A  ■  min(6^, 62) . 

Lemma  2.  L§£  xQ  £  yx  £  y2  £  xx  such  that  r(y1,y2)  and 

yi 

t(yx,y2)  are  inYSEtihi§-  Then  Gy2  is,  a., bilection  (one-one 

correspondence)  with  range  and  domain  X^.  Writing 

yi 

((UpV]),  (u2,v2))  e  Gy*,  we  have  the  formulas 


”uo 

t  -  pt-1  r 

Pt"1 

ur 

(2.3) 

2 

_v2 

■1 

-T-Xr 

t"1 

i 

_vi_ 

and 

(2.4) 


■  i  -1 1 

_  — - 

ui 

rt 

i 

1 

rt 

i 

■o 

u2 

■ 

-1  -1 

_vi. 

rt  t  -rt  P 

7? 

> 

where  r  -  r(y1,y2),  t  -  t(yr,2),  etc. 
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The  proof  follows  immediately  from  Definition  1. 

Y 

Theorem  1.  G  is  a  bi lection  with  range  and  domain 

,  2  5,2 
'£i  is r  all  (y-^yj)  e  D. 

Proof .  Let  (y^,  y2>  be  given.  Partition  the  interval 

[yv  y2J  with  a  8et  where  y1  -  \  <  r]l<  ...  \  -  y2, 

so  that  <  A  for  1  £  i  £  n,  with  A  picked  as 

in  Lemma  1.  Thus  t ^l.)  aiu*  are  Avertible 

^1— 1 

for  1  £  i  £  n.  By  Lemma  2,  is  a  bijection  for  1  £  i  £  n, 

\  ^i.l  %  1  % 

Since  G-,  *  G-,  G_  ,  one  establishes  that  G-  is  a 


yl  11 

1  m  fi 


bijection  for  all  i,  and  in  particular  Gy1  *  G^0  is  a 
bijection. 

Theorem  1  states  that  if  we  are  given  the  fluxes  in 
both  directions  at  a  point  y  in  a  rod  in  which  a  transport 
process  is  taking  place ,  then  we  can  uniquely  determine  the 
fluxes  at  any  other  point  in  the  rod.  We  now  state  the  fact 

that  in  any  transport  process  with  given  incident  fluxes. 

\ 

there  exist  unique  fluxes  at  each  interior  point  of  the  rod. 

U21  |U{| 

Corollary .  Let  -  SCy^yj)  y*  ;  then  for 
each  y  e  Iypy2J  there  exist  unique  vectors  u(y)  e  X 
and  v(y)  €  X  such  that 


- » - -  - - - - - - . . . 
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Proof.  Existence  follows  from  axiom  (iii),  and 
uniqueness  follows  from  Theorem  1. 

Theorem  2.  T(yi>y?)  i™*  c(ypy2)  are  invertible  for 

ail  (ypy2)  €  J  . 

Proof.  Let  (ypy2)  be  given,  and  pick  an  arbitrary 

2 

od  €  X  •  Since  (0,  cp)  e  X  ,  by  Theorem  1  there  exists  a 

V 

unique  (u,  v)  such  that  ((0,  cp) ,  (u,  v))  €  G  .  But  this 
means  that  cp  *  T(y^,y2>  v,  which  establishes  that  T(y^,y2> 
is  invertible. 

Corollary.  The  equations  (2.3)  and  (2.4)  are  valid 
for  all  (yp  y2)  €  D. 

3.  THE  ALGEBRAIC  FORMALISM 

The  family  of  operators  ^ 

*  ft-1  -t-1p 

(3.D  S(yi,y2)  -  -1  T_rt-1„  - 

where  r  *  r(y^,y2),  etc.,  was  defined  by  Redheffer  [5]  in 
order  to  linearize  the  relationship  between  the  scattering 

A 

operators.  The  sigi.  .cance  of  S  is  obtained  from  equation 
(2.4),  which  becomes 


"ii  .  rv 

vJ  ■  s(yi-y2>  [j2J  ’ 
wheie  u^  and  v^(i  *  1,2)  satisfy  the  equation 


*  S(ypy2) 
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Furthermore,  it  is  clear  from  (2—3)  that  S(y^,y2)  Is 
invertible  (for  any  (ypYj'  €  D),  and  its  inverse  is  given 


t  -  Pt 


(3.2) 


-T^r 


Additionally,  we  can  consider  the  "hat"  operator  abstractly 
by  writing 


A  B 


C  D 


-A  1B 


D-CA“1B 


which  is  defined  whenever  A  is  invertible.  It  is  then  easy 
to  verify  that 

A 

A 

s  -  s. 

A—1 

Defining  $  ■  S  ,  we  similarly  have 
5  -  S. 

We  can  npw  express  axiom  (iii)  in  terms  of  the  composition 
of  linear  operators.  From  equation  (2.4),  we  obtain  the 
functional  equations 

A  A  * 

(3.3)  S(y1,y3)  -  S(ypy?)  S(y2,y3)  , 

^yl'y3>  ■  3(y2iy3)  S(yi,y2)  * 

foi  x0  ^  yl  *  y2  *  y3  *  X1  * 
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Theorem  3.  Let  Xq  £  £  y£  £  y$  £  xp  and  write 

Sa  "  S^yl,y2^'  %  "  S^y2,y3^  *  and  Sap  "  S^yl,y3^  * 

Then  I  —  r„P  and  I  -  p  r„  are  Invertible,  and 
-  P  a  -  a  p  - f  . 


<W  • 
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Equations  (3.4)  are  the  algebraic  equations  discovered 
by  Redheffer  [5].  We  can  now  express  our  functional  equation 
in  the  form 

s(yry3)  -  S(y1,y2)  *  S(y2,y3), 
where  the  *— multiplication  is  defined  by 


P'  +  t'P(I  -  r'P)  1  t ' 
t (I  -  r'P)”1  t' 


We  remark  that  the  *— multiplication  is  associative,  and 
that 


so  that  is  a  two-sided  identity.  The  space  of 

fey  2 

all  bounded  linear  operators  onJC  together  with  the  *— 


operation  forms  a  local  semigroup  with  identity.  The 
product 


is  defined  whenever  I  -  B^C2  is  invertible.  (We  also 
need  the  invertibility  of  I  -  but  ex^stence  °£ 

one  inverse  implies  the  other,  for  I  +  C9(I  -  is 

easily  verified  to  be  the  inverse  of  1  -  C2B^.) 

4.  THE  DIFFERENTIAL  EQUATIONS 

Let  A  be  a  function  from  a  real  interval  into  the 
space  of  bounded  linear  operators  on  a  Banach  space.  Then 
A  is  defined  to  be  differentiable  at  if 

lim  r[A(§  r  h)x  -  A(§  )x] 
h-0  h  °  ° 


exists  for  all  x  in  the  Banach  space.  The  derivative 
A' (?Q)  Is  defined  by  equating  A' (§,  )x  to  this  limit.  We 
say  that  A  is  differentiable  on  a  set  Q  if  A  is  differentia¬ 
ble  at  every  point  of  0.  (We  are  defining  u  -ferentiability 
in  the  strong  operator  topology.) 

We  now  state  without  proof  the  following  fact  about 
bounded  linear  operators; 

Theorem  4.  Let  A  and  B  be  functions  mapping  a  real 
interval  into  the  space  of  bounded  linear  operators  on  a 
Banach  space.  If  A  and  B  are  continuous  (in  the  uniform 
operator  topology)  on  the  interval  and  differentiable  (in 
the  strong  operator  topology)  on  a  subset  Cl  of  the  interval, 
we  can  then  state  the  following: 

(i)  AB  is  continuous  (in  the  uniform  operator 
topology)  and  differentiable  on  n  (in  the 
strong  operator  topology) t  and  the  derivative 
of  AB  is  given  by  A*B  +  AB1. 
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(ii)  If  A(?)  Is  also  Invertible  for  all  ?  in  the 
interval,  then  A ~  is  continuous  (in  the 
uniform  operator  topology)  and  differentiable 
on  n  (in  the  strong  operator  topology),  and^ 
the  derivative  of  A-  is  given  by  -  A""  A'A-^. 

Theorem  5.  Suppose  that  the  (one-sided)  derivatives 

y-y  +  h>lh-o  |j;  i 

a  given  interval  and  that 


S(y,y  +  h)|h-Q  and  ^  S(y-h,y)|hiB0  _exist_for  all  y_in 


Ik  +  h> I h-o  *  Ik  s<y  -  h^y>lh-o 


b1(y) 

c(y) 


a(y) 

b2(y) 


Then,  for  fixed  yQ,  the  function  S(y)  -  S(yQ,y)  lg  djf ferm¬ 
ent  iable  on  the  given  interval  and  satisfies  the  differential 
equation 


(4.1) 


S'  - 


+  Pc)  t  a+b]P  +  pb2  +  Pep 

ket  T(b2  +  cP) 

t  p" 
r  t_ 

The  proof  follows  immediately  from  the  preceding  theorem 
and  the  functional  equations 


where  we  let  S  * 


as  usual. 


and 


S(y0,y  +  h)  -  S(yQ,y)  *  S(y,  y  +  h) 

S(y0,  y  -  h)  -  S(yQ,y)  3(y  -  h,y), 

* 

3  (y0,  y  -  h)  -  S(y  -  h,y)  S(yQJy)  . 

Equation  (4.1)  is  the  well-known  invariant- imbedding 
equation  (see  [5 J)  and  is  identical  to  equation  (1.2). 

The  equations  (1.1)  for  the  "internal  fluxes"  also  follow 
Immediately  from  the  equations 


—>20*” 


u(y  +  h) 

u(yf| 

v(y  +  h) 

j  *  S(y,  y  +  *0 

_v(y)J 

u(y  -  h)| 

A 

u(yf 

v(y  ~  h) 
— 

-  §(y  -  h,y) 

v(y) 
L»  J 

s.  application 

We  now  indicate  an  application  of  the  operator 
differential  equation  (1.2).  (Other  applications  can  be 
found  in  [1,6]).  We  consider  a  one-dimensional  neutron- 
transport  process  in  which  the  energy  of  the  neutrons  is 
allowed  to  vary  continuously.  The  following  functions 
shall  represent  the  physical  parameters  that  describe  the 
neutron  transport  process: 


a^(y,  u)|  the  reciprocal  of  the  mean  free  path  o 

a  (y,  u)[  a  neutron  of  energy  U  travelling  to  the 

"1  at  the  point  y  in  the  rod. 


jrightj 


f^(y,  v>| 

fr(y>  v>| 


g^(y>  v, 


the  expected  energy  distribution  (in  u) 
of  neutrons  travelling  to  the  jrightj 
due  to  a  fission  at  y  of  a  neutron  or 

u  Jleft  1 

energy  v  travelling  to  the  * 

the  expected  energy  distribution  (in  u) 

,  fright 

of  neutrons  travelling  to  the  |ieft  J 

due  to  a  fission  at  y  of  a  neutron  of 

.  fleft  \ 

energy  v  travelling  to  the  ’‘jjightj 
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If  we  represent  neutron  fluxes  as  energy-distribution 
functions  u(v)  and  v(v;>  the  above  definitions  yield  the 
integro'-differential  equations 

|y  u(y,  n)  -  J  or(y,  v)fr(y,  \i,  v)u(y,  v)dv 

-v  (y,  n)  u(y,  u) 

+  J  °i(y>  v)g^(y,  v)v(y,  v)dv, 

-  v(y,  n)  -  J  or(y,  v)gr(y,  n,  v>u(y,  v)dv 
+  J  o^(y,  v)f^(y,  u,  v)v(y,  v)dv 

-  o^(y,  u)v(y,  n), 

where  the  integrations  are  over  the  allowable  energy 
range . 

Comparing  equation  (5.1)  with  (1.1),  we  see  that  the 
operators  a(y),  b^(y),  b2<y),  c(y)  are  given  by 

a(y)  y  (n)  -  J  ot(y,  v)g^(y,n,v)  Y(v)dv, 
c(y)  ¥  (u)  -  J  or(y,  v)gr(y,^,v)  ¥(v)dv  , 

(5.2) 

bx(y)  ¥(u)  «  J  ar( y,  v)fr(y,^,v)  ¥0>)dv 

-  or(y,u)  *  00, 

b2(y)  ¥(u)  -  J  at(y,v)f^(y,n,v)  ¥(v)dv 

-  ^(y,^)  *00* 

where  ¥  is  an  energy  distribution  function  representing  a 
neutron  flux. 
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We  then  represent  our  transport  operator'  as  follows: 


P(x)  *  (u)  -  J*  R(x,  n,  v)  f  (v)dv, 

r(x)  ijr  (u)  ■  J  r(x,  u,  v)  i|i  (v)dv, 

(5.3) 

t(x)  i|i  (u)  -  Tq(x,  |i)  *(u)  +  J*  T(x,  [x,  v)  *(v)dv, 

t(x)  ii  (^)  -  tQ(x,  n)  ijr(u)  +J  t(x,  m,  v)  *(v)dv. 

We  can  now  apply  our  operator  differential  equations 
(1.2)  to  obtain  the  differential  equations  for  the  "reflec¬ 
tion"  and  "transmission"  functions  defined  by  (5.3);  thus: 


V) 


°l(x>  v)g^(x,  U,  v) 


(5.4) 


h  Vx> 

k  to<x<  »> 
T(x,  u,  v) 


~  [ot(x,  v)  +ar(x,  n)]R(x,  n,  v) 

+  J*  v)f^(x,  a,  v)R(x,  a)da 

+  J  or(x,  “)fr(x,  U,  a)R(x,  a,  v)da 
+  JJ*  R(x,  u,  3)or(x,  a)gr(x,  3,  a)R(x,  a,  v)dad3, 
-  -  0^(x,  u)To(x,  n)  , 

■  “  <*r(x,  n)tQ(x,  u)  , 

"  J  T(x>  a)o^(x,  v)f^(x,  a,  v)da 
""  T(x,  n,  v)a^(x,  v) 

+  n  t(x,  a)°T(x,  P)gj.(x,  a,  3)R(x,  3,  V)dad3' 
+  Tq(x,  \i)[ot(x,  v)f^(x,  »,  V)  I- J  ar(x,  a) 
gy  »Xj  *i,  a)R(x,  a,  v)da]  , 
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|j  t(x,  Ut  v)  ■  J  or(x,  a)fr(x,  u,  “)t(x,  a,v)da 

-  ar(x,  n)t(x,  v) 

+  JJ  R(x,  n,  a)or(x,  P)gr(x,  a,  P)t(x,  P,  v)dadp 

+  t  (x,  v)  [a  (x,  v ) f  (x ,  a,  v) 
o  r  r 

+  J  R(x,  n,  a)ar(x,  vjg^x,  a,  v)da]  , 

|-  r(x,  u,  v)  »  JJ  T(x,  n,  a)ar(x,  ^)gj.(x,  a,  P)t(?  ,  P,  v)dadp 

+  Tq(x,  n)  J  ar(x,  aOg^x,  n,  a)t(x,  a,  v)da 
+  tQ(x,  v)  J  T(x,  n,  a)ar(x,  v)gJ.(xJ  a,  v)da 
+  \i,  v)TQ(x,  u)tQ(x,  v). 

The  Initial  conditions  are 

R(0,  \i,  v)  -  r(0,  |i,  v)  -  T(0,  n,  v)  «  t(0,  n,  v)  -  0, 

(5.5)  V  , 

TQ(0,  n)  -  tQ(0,  \x)  -  1. 

We  see  that  TQ(x,  n)  gives  the  fraction  of  the  neutrons 
of  energy  \i  incident  to  the  left  that  leave  the  reactor  withr- 
out  entering  into  any  fission  reactions,  and  tQ(x,  n)  rep¬ 
resents  the  similar  fraction  for  neutrons  incident  to  the 
right.  While  these  rather  formidable-looking  equations  are 
difficult  to  obtain  directly  from  physical  considerations, 
they  follow  from  our  operator  equations  (1.2)  with  sur¬ 
prisingly  little  mathematical  manipulation- 
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6a _ BEMABKS 

We  have  given  an  axiomatized  description  of  one¬ 
dimensional  transport  processes  and  have  shown  how  our 
axiomatized  system  leads  to  the  desired  algebraic  and 
differential  equations  of  transport  theory.  One  could 
also  start  with  the  differential  equations  (1.2)  of  the 
transport  process  and  then  verify  the  algebraic  equations 
(3.4).  This  approach  has  been  used  by  Redheffer  [8]. 

In  Sec.  5  we  used  the  generalized  invariant-imbedding 
equctions  (1.2)  to  describe  an  energy-dependent  neutron- 
transport  process.  We  remind  the  reader  that  when  we 
consider  neutron  transport  processes,  r \xr  formalism  deals 
only  with  rods  of  less  than  "critical  length."  When  a 
rod  reaches  critical  length,  an  infinite  number  of  neutrons 
is  produced,  and  the  fluxes  can  no  longer  be  considered 
as  elements  of  the  Banach  space  X  •  To  compute  the 
critical  length  by  the  invariant— imbedding  method,  we 
find  the  point  where  the  solution  of  the  first  equation 
of  (1.2)  becomes  "infinite,"  as  described  in  [1]  and  [3]. 

We  conclude  by  remarking  that  our  mathematical  method 
of  approach  in  this  Memorandum  is  based  on  the  idea  of 
applying  the  invariant— imbedding  concept  to  fand  1 i of 
operators .  It  seems  likely  that  this  technique  may  have 
further  applications  in  the  domain  of  mathematical  physics. 
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